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^ . Abstract: 



The complete electroweak 0(a) radiative corrections to the Higgs-boson production 
processes e + e~ — > V[ViH (I = e, /i, r) are calculated in the electroweak Standard Model. 
For e + e~ — > z/ e ^ e H, where ZH production and W-boson fusion contribute, both production 
channels are added coherently. The calculation of the corrections is described in some de- 
tail including, in particular, the treatment of the Z-boson resonance in the ZH-production 
channel. The discussion of numerical results focusses on the total cross section as well as 
on angular and energy distributions of the Higgs boson. In the G^-scheme, the bulk of 
the corrections is due to initial-state radiation. The corrections turn out to reduce the 
total cross section by ~ 10% for high energies, where the W-boson fusion dominates. In 
this region, the corrections depend only weakly on the energy and the production angle 
of the Higgs boson. Based on an analysis of the leading universal corrections, a simple 
improved Born approximation is introduced. This approximation describes the corrected 
cross section within about 3%. 
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1 Introduction 



One of the most important open problems of particle physics is the understanding of 
the mechanism of electroweak symmetry breaking. In the electroweak Standard Model 
(SM) it is provided by the Higgs mechanism, leading to the prediction of a physical scalar 
particle, the Higgs boson. The investigation of the mechanism of electroweak symmetry 
breaking and, in particular, of the Higgs boson, will be one of the main concerns at the 
Large Hadron Collider (LHC) at CERN. The LHC experiments ATLAS [ d and CMS 
[ |2] are sensitive to the SM Higgs boson over the whole mass range from the present 
lower experimental limit of 114.4 GeV [0] up to 1 TeV and will discover the Higgs boson, 
if it exists and has no particularly exotic properties. Moreover, these experiments will 
determine various properties of the Higgs boson, such as its mass, branching ratios, and 
ratios of its coupling constants. 

However, the complete profile of the Higgs boson can only be studied in the clean 
environment of an electron-positron linear collider [ EJ E] • In e + e~ annihilation there 
are two main production mechanisms for the SM Higgs boson. In the Higgs-strahlung 
process, e + e~ — > ZH, a virtual Z boson decays into a Z boson and a Higgs boson. The 
corresponding cross section rises sharply at threshold to a maximum of a few tens of GeV 
above Mi + Mh and then falls off as where \fs is the centre-of-mass (CM) energy 
of the e + e~ system. In the W-boson-fusion process, e + e~ — > i^eH, the incoming e + and 
e~ each emit a virtual W boson which fuse into a Higgs boson. The cross section of the 
W-boson-fusion process grows as In s and thus is the dominant production mechanism for 
y/s ^> Mh. The cross section for the similar Z-boson-fusion process, e + e~ — > e + e~H, is 
about one order of magnitude smaller. 

At a linear e + e~ collider with a CM energy of 500 GeV and an integrated luminosity 
of 500 fb~ , of the order of 10 4 Higgs bosons can be produced per year [E]. This allows to 
measure the Higgs-production cross sections and thus the Higgs-gauge-boson couplings 
at the level of a few per cent. Consequently, adequate theoretical predictions have to 
take into account radiative corrections and the relevant effects of the finite decay width 
of the Z and Higgs bosons. For a heavy Higgs, which decays mainly into W-boson pairs, 
finite-width effects have been investigated in Ref. [ISJE]- On the other hand, the effects 
of the finite width of the Higgs boson can be neglected if the Higgs boson is light, i.e. if 
its width is small. 

The Higgs-strahlung process has been investigated in lowest order in Ref. [I10j. In- 
cluding the Z-boson decay, the process e + e~ — > Z*H — > //H has been studied by Bjorken 
[ lllj . For the process e + e~ — > ZH, the 0(a) electroweak corrections have been calculated 
in the soft-photon approximation in Refs. [ EH , and a Monte Carlo algorithm for 
the calculation of the real photonic corrections to this process was described in Ref. [Ho]. 
A compact analytical formula for the electromagnetic corrections to the total cross section 
can be found in Ref. [I13j. 

The vector-boson- fusion processes have been investigated in lowest order in Ref. [ . 
The electroweak corrections to e + e~ — > vvH have attracted a lot of interest recently. Ana- 
lytical results for the one- loop corrections to this process have been obtained in Ref. [ E| 
as MAPLE output, but a numerical evaluation of these results is not yet available. A 
first complete calculation of the 0(a) electroweak corrections to e + e~ — > z/z/H in the SM 
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has been performed in Ref. [QE]- The contributions of fermion and sfermion loops in 
the Minimal Supersymmetric Standard Model have been evaluated in Refs. [ EH 120] with 
seemingly differing results. We have performed a completely independent calculation of 
the 0(a) electroweak corrections to the complete process e + e~ — > vvH in the SM. First 
results of our calculation have already been presented in Ref. [|2T]. There, we have suc- 
cessfully compared our results for the total cross section with those of Refs. [ [TBI EH 
and pointed out that the differences between Refs. [ ITU] and [ 1201 are due to different 
renormalization schemes and input parameters. 

In this paper, we present the details of our calculation of the 0(a) electroweak cor- 
rections to the processes e + e~ — > z/ e z/ e H, z/^^H, and i/ T z/ T H. While the first process gets 
contributions from W fusion and Higgs-strahlung, the final states with \x or r neutrinos 
receive contributions only from Higgs-strahlung off Z bosons. Single hard-photon radia- 
tion is included using the complete matrix element, and higher-order ISR corrections are 
taken into account in the leading-logarithmic approximation. The finite Z-boson decay 
width is introduced in the constant-width scheme and in a gauge-invariant scheme that is 
based on a factorization of the Z resonance in the gauge-invariant set of diagrams related 
to the (neutral-current) Higgs-strahlung process. 

The paper is organized as follows: in Section |2] the calculation of the virtual, real, 
and (higher-order) leading-logarithmic corrections is described, and the leading universal 
corrections are discussed. Section El contains a discussion of numerical results. The paper 
is summarized in Section HJ Further useful information on the calculation is collected in 
the Appendix. 

2 Calculation of radiative corrections 

2.1 Conventions and lowest-order cross section 

We consider the processes 

e~(px,ai) + e + (p 2 ,a 2 ) — ► ui(k{) + i>i(k 2 ) + H(fc 3 ), l = e,fi,r, (2.1) 

where the momenta Pi, kj of all particles and the electron helicities <7j are given in paren- 
theses. The helicities take the values oi = ±1/2, but we often use only the sign to 
indicate the helicity. The electron mass is neglected whenever possible, i.e. it is kept 
finite only in the mass-singular logarithms related to initial-state radiation. This im- 
plies that the lowest-order and one- loop amplitudes vanish unless <j\ = —a 2 . There- 
fore, we define a — a± — —a 2 . The particle momenta obey the mass-shell conditions 
Pi = v\ = = k 2 = an d k\ = My. For later use, the following set of kinematical 
invariants is defined: 

8 = (pi+p 2 ) 2 , 

Sij — {k% + kj) , i,j = 1, 2, 3, 

Uj = (Pi - kj) 2 , % = 1, 2, j = 1, 2, 3. (2.2) 

In lowest order the processes (|2.1j) proceed via the diagrams shown in Figure More 
precisely, only for electron neutrinos in the final state (/ = e) both the ZH-production and 
WW-fusion diagrams contribute, while for /i and r neutrinos merely the former exists. 
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Figure 1: Lowest-order diagrams for e e + — > uuH 



In the calculation of the tree-level amplitude Ml and of the one-loop amplitude Ml, 
which is described in the next section, we separate the fermion spinor chains by defining 
standard matrix elements (SME). To introduce a compact notation for the SME, the 
tensors 

^{Zafh} = ^+(P2) {7a,7a7/37 7 }^ M e-(Pl), 
F {a,a/3 7 } = U Vl {kl) {7a,7a7/37 7 } U-V Pl (k 2 ), 
r {a,a/3 7 } = U n (h) {7 Q ,7 Q 7 /3 7 7 }w_M c -(pi), 

T {a,af3y} = (p 2 ) {7a, 7a7/37 7 } ^-^ (2.3) 

are defined with obvious notations for the Dirac spinors v c +(p 2 ), etc., and uj± = (1 ±7 5 )/2 
denote the right- and left-handed chirality projectors. Here and in the following, each 
entry in the set within curly brackets refers to a single object, i.e. from the first line in 
the equation above we have r^ e,fT = v c +{p2)la^ 1 aU c -{pi) , etc. Furthermore, symbols like 
T p are used as shorthand for the contraction T^p^. For the ZH-production channel we 
define the 26 SME 

A/fZH,(T pee, a T^vu,{ot,apip2} A/f^ ,cr pcc,<r ~puu,{a,apip2} 

J '{1,2} ~~ 1 a 1 ? 7 {3,4} ~~ 1 afcifc 2 1 ' 

a^ZH,(t _ -pec,o- -n^i/,{pi,p 2 } /Ci ZR ' CT — F ee ' cr F^lPi^} 

- /Vt {5,6} — 1 fci 1 5 y l {7,8} — 1 fc 2 1 ' 

«"/ZH,(j _ pee, o- -puu,{a/3pi,aj3p2} A/f 211 '"" — F ee,fJ Y vv A a ftVli a ftVl\ 

yvl {9,10} ~~ 1 1 ) • /VI {11,12} ~~ 1 a/3k 2 1 ' 

yVl™'' 7 = r w -*. (2.4) 

For the WW-fusion channel we introduce the following set of 13 SME, 

A/tfWW ~pve peiA{a,afci»i} A/fWW ~pue pe^,{o,afcipi} 

^{l^} - L a L i yW {3,4} - i afc 2 P2 1 ' 

UWW _ pi/e pei^,{pi,fci} A/f ww — F^ e F el/ '{Pi> fc i} 

yvl {5.6} ~~ P2 1 ' y l {7,8} — i fc 2 1 ' 

A^fWW -pue -pev,taf3pi,af3ki} uWW ~Tve -peu,(a/3pi,afSki} 

JVl {9,10} ~ 1 af3p 2 1 ' - /Vt {ll,12} — i a/3fe 2 1 ' 

/VfWW = pei/,007. (2.5) 

The tree-level and one-loop amplitudes can be expanded in terms of linear combina- 
tions of SME, 

Ml = £ + 5,- £ F^r, n = 0, 1, (2.6) 

i=i i=i 
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with Lorentz-invariant functions F^J ,cr and i^ w - This decomposition is unique in D 
dimensions, i.e. if only the Dirac equation for the spinors is used. The number of SME 
can, however, be reduced further by exploiting the four-dimensionality of space-time, 
which implies relations among the SME given above. In fact, it is possible to express 
the set of all 39 SME in terms of two independent SME. We list the expressions of the 
two independent SME, which can be identified with M 1 ' , and the relations among the 
SME in the Appendix. 

The lowest-order amplitude reads 



where 



Mq = .M ZH ' CT + 5*-M™ w , (2.7) 



M ° ~ " 244 (s - Mf )( Sl2 - M| + iM z r z ) ^ ' (2 ' 8) 



yww _ r e 3 M w 1 A>WW /o q\ 



with the chiral couplings g% c of the electron to the Z boson, 



The sine and cosine of the weak mixing angle are fixed by 

4 = i-4 = |f- (2.H) 

For s 12 — > M|, the lowest-order amplitude develops a resonance corresponding to real ZH 
production with the subsequent Z — > z/^ decay. We, therefore, have included a finite 
Z-boson width T z in the denominator, which results from a (partial) Dyson summation of 
the corresponding propagator. This procedure is discussed in more detail in the context 
of the virtual radiative corrections in the next section. 
Finally, the lowest-order cross section reads 

ao = hl m £ I^ 1 + * P -°)Q - 2P + a ) l-^ol 2 , (2-12) 

where P± are the degrees of polarization of the e* 1 beams and the phase-space integral is 
defined by 

/ d *> = (A/tlw) (^(p 1+ P 2 -g%). (2.i3) 

2.2 Virtual corrections 

2.2. 1 Survey of one-loop diagrams 

The virtual corrections can be classified into self-energy, vertex, box, and pentagon 
corrections. The generic contributions of the different vertex functions are shown in 
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Figure 4: Diagrams for the z/^ZH and v\Vffift. vertex functions 

Figure El The first three lines contain those diagrams that contribute to all z/^H final 
states, whereas the diagrams in the last three lines contribute only to e + e~ — > i/ e i> e H. 

The complete set of pentagon diagrams is shown in Figure H3 The last eight diagrams 
contribute only for the z/ e z/ e H final state. The z/^ZH and box diagrams are depicted 

in Figure El the e + e~ZH box diagrams in Figure and the e~u e W + B. box diagrams in 
Figure |H1 Those for the e + f e W _ H box diagrams can be obtained by charge conjugation. 

The diagrams for the i^/H and e + e~H vertex functions are listed in Figure Figure[S] 
shows the Feynman diagrams for the W + W~H vertex function and Figure 01 those for the 
ZZH and 7ZH vertex functions. Note that in Figure 01 those diagrams that are obtained 
from the diagrams in the first three lines of this figure by reversing the charge flow in 
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Figure 5: Diagrams for the e + e ZH vertex function 
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Figure 9: Diagrams for the ZZH and 7ZH vertex functions 
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the loop have been suppressed. Most of the diagrams for the self-energies and the Vxv{L, 

e + e~Z, and e ± fjw vertex functions can be found in Ref. [122]. 

All pentagon and box diagrams are UV finite, and also the e + e~H and viViH. vertex 
functions are finite since we neglect the electron mass everywhere apart from the mass- 
singular logarithms. For the other vertex functions, 7ZH, ZZH, Vip{L, e + e~Z, W + W~H, 

e^e W, and for the ZZ, 7Z, and WW self energies the corresponding counterterm dia- 
grams have to be included. 

2.2.2 Calculational framework 

The actual calculation of the one-loop diagrams has been carried out in the 't Hooft- 
Feynman gauge using standard techniques. The Feynman graphs have been generated 
with FeynArts [ 23 j and are evaluated in two completely independent ways, leading to 
two independent computer codes. The results of the two codes are in good numerical 
agreement (i.e. within about 12 digits for non-exceptional phase-space points). Apart from 
the 5-point integrals, in both calculations the tensor coefficients of the one-loop integrals 
are algebraically reduced to scalar integrals with the Passarino-Veltman algorithm [ 12*1] 
at the numerical level. The scalar integrals are evaluated using the methods and results 
of Refs. [ Em EI) , where ultraviolet divergences are regulated dimensionally and IR 
divergences with an infinitesimal photon mass m 7 . The renormalization is carried out in 
the on-shell renormalization scheme, as e.g. described in Ref. [I27j. 

In the first calculation, the Feynman graphs are generated with FeynArts version 1.0 [ 
|2*3*] . With the help of Mathematica routines the amplitudes are expressed in terms of SME 
and coefficients of tensor integrals. The output is processed into a Fortran program for the 
numerical evaluation. For the evaluation of the tensor 5-point function two approaches 
have been followed: the usual Passarino-Veltman reduction and the direct reduction to 
4-point integrals of Ref. [I2H]- The results based on the Passarino-Veltman algorithm 
become numerically unstable at the phase-space boundary and could only be rescued by 
a careful extrapolation out of the numerically safe inner phase-space domains, as described 
in Section 2.2.4 of Ref. [ 29, . The instabilities are due to the occurrence of inverse Gram 
determinants in the recursive tensor reduction. The direct reduction of Ref. [ |2E] avoids 
such inverse Gram determinants, rendering the results of this approach well behaved near 
the phase-space boundary. 

The calculation of the virtual corrections has been repeated using the background- 
field method [EU], where the individual contributions from self-energy, vertex, box, and 
pentagon corrections differ from their counterparts in the conventional formalism. The 
total one-loop corrections of the conventional and of the background-field approach were 
found to be in perfect numerical agreement. 

The second calculation has been made using FeynArts version 3 [ HH] for the gener- 
ation and FormCalc [E2] for the evaluation of the amplitudes. The analytical results of 
FormCalc in terms of SME and their coefficients were translated into C++ code, and 
the interference of the one-loop with the lowest-order amplitude calculated numerically. 
However, in the case of the pentagon diagrams this interference was calculated analyt- 
ically using FeynCalc [ 33J. After evaluation of the fermion traces, the loop momenta 
appear in the numerator only in scalar products and can be cancelled against propagator 
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denominators. In this way, only scalar 5-point functions remain, thereby avoiding inverse 
Gram determinants and the related numerical instabilities. 

Finally, the contribution of the virtual corrections to the cross section is given by 

Sa viTt = i- J d$ 3 J2 \{l + 2P-a){l-2P +( r)2Re{M<i{Miy}. (2.14) 



2.2.3 Treatment of the Z-boson resonance 



At tree level the introduction of the Z-boson decay width in the lowest-order matrix el- 
ement A^o of Eq. (|2.7|) did not pose any problem with gauge invariance or double-counting 
of higher-order effects, since the gauge-invariant ZH-production part before Dyson sum- 



mation, Ai^ ,u (T z = 0), was simply rescaled as 



Sl2 



Ml 



M% n ' a (r z = o). 



Sl2 



- M| + iAf : 



(2.15) 



z- 1 - z 



This modification reproduces the correct Breit-Wigner shape near the resonance and 
changes Mq ,a {T z = 0) at the relative order 0{T z /M z ) ~ 0(a) away from the resonance. 
It should also be mentioned that we introduce a fixed Z-boson width (in contrast to a 
running width), i.e. the Z-boson mass M z deviates from the on-shell mass at the two- loop 
level (see e.g. Ref. [ED- 

The issue of gauge invariance, resonances, Dyson summation, and radiative corrections 
is rather complex. In fact, no simple general solution for a gauge-invariant treatment of 
finite-width effects exists yet. Fortunately, the situation in our case is relatively simple. 
We proceed in two different ways. 

First we apply the so-called naive fixed-width scheme which simply means that each 
resonance propagator l/(si 2 — M|) is replaced by l/(si 2 — M|+iM z r z ), while non-resonant 
contributions are kept untouched. The contribution iM z r z originates from the imaginary 
part of the (transverse part of the) Z self-energy £ zz on resonance; in order to avoid 
double-counting, the 0{a) contribution thus contained in A4 ZH,<T has to be subtracted 
from the one-loop amplitude. In summary, the one-loop amplitude in the fixed-width 
scheme reads 



fixed- width 



Ai ZH ' CT (r z = o)-^ z . self (r z = o) 



f| s 12 -M|+iM z r z 



Ef(s) - Re{S zz (M|)} Ef(s 12 ) - Sf (M| 



s-M| 



+ 



S12 - Mi 



+ 26Z ZZ 



M 



ZH> 




(2.16) 



where we have made the Z self-energy and its corresponding counterterm contributions 
explicit. Note that the term Re{S| z (M|)} is identical with the mass counterterm 6M Z , 
while S| Z (M|) receives its imaginary part iIm{S| z (M|)} in the resonant part from the 
subtraction of the iM z r z contribution already contained in A4q R,(T . Since only complete 
S'-matrix elements exhibit the gauge-invariance properties such as the independence of 



11 



gauge-fixing conditions (gauge parameters) and the validity of Slavnov-Taylor identities, 
this procedure potentially violates gauge invariance. 

As a second option, we applied a factorization scheme where the full ZH-production 
amplitude before Dyson summation is rescaled similar to A4q ' a (Tz = 0) in Eq. (j2.15|) . 
This procedure is analogous to the treatment of the W-boson resonance in pp — ► W — > Ivt 
as described in Ref. [EHj. Of course, also here double-counting of iM z r z terms has to be 
avoided. In summary, the one-loop amplitude in the factorization scheme reads 

- s 12 -Mi + iM z T z M ' (Fz - 0) + s 12 -M 2 M ° 



s X2 -M£ + iM z r z 



£ zz (s) - Re{E zz (M|)} S zz (si 2 ) - S ZZ (M, 



Z/ 



2 + : m ^ 25z ™ 



s-M4 S12 - M 



z 





(2.17) 



Note that this procedure preserves gauge invariance, since the gauge-invariant amplitude 
A / l ZH ' cr (r z = 0) is only rescaled and accompanied by another gauge-invariant term propor- 
tional to A / f ZH '°\ However, the rescaling puts all non-resonant terms in A4 ZH,<T (r z = 0) to 
zero on the resonance. The corresponding error is of the order of the non- resonant terms 
in the 0(a) corrections, i.e. of order C(ar z /M z ) ~ 0(a 2 ). 

Within integration errors, both schemes give the same results. 

2.2-4 Universal electroweak corrections 

The electroweak corrections contain large contributions of universal origin. Besides 
initial-state radiation (ISR), which is discussed in Section \2 .3 .3| these consist, in partic- 
ular, of the corrections associated with the running of a and corrections proportional to 
m 2 /M^j that can be associated to the p parameter or to the renormalization of the elec- 
troweak mixing angle. By suitable parametrization of the lowest-order matrix elements, 
these universal corrections can be incorporated into the lowest order thus reducing the 
remaining corrections. This does not only reduce the 0(a) corrections but in general also 
the higher-order corrections. 

The running of u(Q 2 ) from Q 2 = to Q 2 = M| is of the order of Aa ~ 6%. Since 
the cross section for e + e~ — ► vvH is proportional to a 3 this amounts to an effect of 
~ 18% if a(0) defines the electromagnetic coupling ("a(0)-scheme"). By parametrizing 
the lowest-order cross section in terms of a(M|) ("a(M|)-scheme") this large effect can be 
incorporated into the leading-order expressions. Note that Aa contains large contributions 
from the hadronic vacuum polarization that cannot be calculated in perturbation theory. 
Thus, the sensitivity to these effects is absorbed in the lowest-order cross section via 
a(M|) in the a(M|)-scheme. Alternatively the electromagnetic coupling can be deduced 
from the Fermi constant ( "G^-scheme" ) via 

7T 
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where G M is measured in muon decay. Consequently, in the transition from the a(0)- 
to the G M -scheme the constant 3Ar has to be subtracted from the relative correction to 
a cross section that is proportional to a 3 , where Ar contains the electroweak radiative 
corrections to muon decay. Since Ar is about 3% for the empirical value of m t , this 
shifts the relative corrections by ~ 9% with respect to the a(0)-scheme. Apart from 
absorbing Aa, the quantity Ar additionally contains universal corrections proportional 
to ml/My; that can be associated to the p parameter, Ar — Aa ~ — Apc^/s^ ~ —3% 
with Ap = 3G M m?/(8 V2n 2 ) [EE]. 

The cross section for e + e~ — > vvH is dominated by the WW-fusion diagram, which 
gets its main contribution from the region of small momentum transfers. Consequently, 
the corresponding corrections are determined by the eu e W and WWH vertex corrections 
for small invariant W masses and depend only weakly on the energy. The correction to 
the ez/ e W vertex in the relevant kinematical region is well approximated by Ar. It turns 
out that this is also the case for the main contributions to the WWH vertex. Thus, 
parametrizing the lowest order in terms of (G^-scheme) absorbs a large part of the 
universal corrections. Since in the G^-scheme all large universal corrections related to the 
running of a and most of the corrections oc m%/Myj are absorbed, we use this scheme in 
the following. With respect to this scheme, the corrections in the a(0)-scheme are shifted 
by 3Ar « +9% and those in the a(Mf )-scheme by 3(Ar - Aa (M|)) « -9%. 

We have extracted the leading m t -dependent corrections in the heavy-top limit in the 
G,jr scheme: 



ZH,, 



a: 



G M -scheme m t ^oo ^ ng 2 



Myj 3c w s w M w 



M% u,tr 
9z c 



M? w ~ + « (2-19) 

The leading m^-enhanced terms of the WW channel agree with the terms derived for the 
HWW vertex [ EZj , since in the G M -scheme all leading contributions related to the 
W-boson coupling to fermions are absorbed in ckg m - In the ZH channel, m^-enhanced 
terms do not only result from the HZZ vertex, but there are also remnants originating 
from the renormalization of the Z-boson couplings to fermions. In contrast to the WW 
channel, in the ZH channel there are also logarithmic terms lnm t for a large top-quark 
mass. 

In the physically interesting region of e + e~ — > z/z/H, ^/s is of the same order as m t 
or larger. Nevertheless, the expressions for m t — > oo reproduce the full m t -dependent 
corrections rather well for the WW channel, which is dominated by small momentum 
transfers. The situation for the ZH channel, where ^/s is a typical scale for the momentum 
transfer, is completely different. There, the expressions for m t — > oo do not provide a 
good approximation, and we did not succeed in finding a simple approximation for the 
fermionic corrections to the ZH channel. This is due to the presence of loop integrals 
that depend both on s and m%. Since, however, the cross section for e + e~ — ► z/z/H is 
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Figure 10: Feynman diagrams for e + e — > z/z/IFy 

dominated by the WW channel we consider it useful to define the following improved 
Born approximation 

i non-photonic i J^WW ^ a m t ( n on \ 

da IBA =da -da _ ^. (2.20) 

This cross section is then convoluted with the ISR as given in ()2.49|) below. 

As discussed in Ref. [EI] and in Section the bosonic corrections are small for the 
WW channel but large for the ZH channel. In order to find the source of these large 
corrections, we have evaluated the leading bosonic corrections in the limit s ~ 2> 
which are of the order 0(s/M^) ~ 0(M^/M^). We found that these terms are small and 
cannot explain the large bosonic corrections in the ZH channel. In the 't Hooft-Feynman 
gauge the large corrections arise predominantly from box diagrams involving W-boson 
exchange [H3j. The enhancement of these diagrams is partially due to the gauge-boson 
coupling to the electron which for W bosons is stronger than for Z bosons. 

2.3 Real photonic corrections 

2.3.1 Matrix- element calculation 

The real photonic corrections are induced by the process 

e-( Ph ax) +e + G» 2 ,cr 2 ) — > u^kx) + Pi(k 2 ) + H(fc 3 ) + 7^, A), l = e,n,r, (2.21) 

where k and A denote the photon momentum and helicity, respectively. The Feynman 
diagrams of this process are shown in Figure [TU1 We have evaluated the helicity ma- 
trix elements A4° x of this process using the Weyl-van der Waerden spinor technique as 
formulated in Ref. [OH]. The amplitudes for the helicity channels with o~\ = o 2 vanish 
for massless electrons. The finite-mass effects of electrons and positrons are included in 
the treatment of soft and collinear singularities, as described in the next section. In the 
notation of Ref. [EH] the amplitudes A4° x read 
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7 44 (Pi + - A;) 2 - M| (ki + k 2 y - Mi + iM z r z 

x A^( Pl ,p 2 ,h,k 2 ,k), 

WWiA= y/2e 4 M w 1 1 

' e 4 (Pi-fci) 2 -^ (P2-^ 2 ) 2 -M^ 

xirbbfeii^^) (2.22) 
with the auxiliary functions 

^4 (pi,p2,h,k 2 ,k) = -( Pl k 2 ) fo^fofc). > 

A™ x (p uP2 ,k u k 2 ,k) = -A^(p 2 , Pl ,k 2 ,k 1 ,ky, 

A™_ x ( PuP2 , k u k 2 , k) = A™(p 1 ,p 2 , k 2 , h, k)*, 

4 WW/ „ r. r. M /„L\J (PlP2)*(P2h) (fcifc) ((Pifci)*(p2fcl) + (Plk)*( P2 k)) 

A ~ ~ ^ \- { p lk Y{p 2 kY + (pik)* [(pi - fc x - fc) 2 - M w ] 

(fcifc) (fc 2 fc) {(p 2 k 2 y(p 2 h) - (hkyjhk)) \ 

(p 2 k)* (p 2 k)* [(p 2 - fc 2 - fc) 2 - M w ] /' 
^ W A W (Pi,P2,fci,fc 2 ,fc) = -Ar v (p 2 ,p 1 ,fc 2 ,A: 1 ,fc)*. (2.23) 

The relations between the A;;; functions that differ only in the photon helicity A are a 
consequence of the CP symmetry of the process, while the relation associated with a 
reversion of both a and A results from a P transformation. The spinor products (...) are 
defined by 



(pq) = e AB p A qB = Zy/ptffo 



Op . 6q _ i(j> Oq . p 

— sin e Vq cos — sin — 

2 2 2 2 



(2.24) 



where pa, Qa are the associated momentum spinors for the light-like momenta 



P^ — Po(l; si 11 Op cos (f) p , sin P sin p , cos p ) , 

(f — qo{l,smO g cos(f> q ,smO q sm(f) q ,cosO q ). (2.25) 

The matrix elements (|2.22|) . (J2.23)) have been successfully checked against the result 
obtained with the package Madgraph [EH] numerically. 

The contribution cr 7 of the radiative process to the cross section is given by 

a 7 = Wd$ 7 £ hl + 2P_a)(l-2P + a) £ \M° X \ 2 , (2.26) 

/S cr=±i A=±l 

° 31 2 

where the phase-space integral is defined by 

/ d * W p|b (n / (» + - * - g *») • < 2 - 27 > 
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2.3.2 Treatment of soft and collinear singularities 



Without soft and collinear regulators the phase-space integral (|2.2fij) diverges in the 
soft (fco — > 0) and collinear (jpik — > 0) phase-space regions. In the following we describe two 
procedures of treating soft and collinear photon emission: one is based on a subtraction 
method, the other on phase-space slicing. In both cases soft and collinear singularities 
are regularized by an infinitesimal photon mass and a small electron mass, respectively. 

(i) The dipole subtraction approach 

The idea of so-called subtraction methods is to subtract a simple auxiliary function 
from the singular integrand of the bremsstrahlung integral and to add this contribution 
back again after partial analytic integration. This auxiliary function, denoted |-M S ub| 2 
in the following, has to be chosen in such a way that it cancels all singularities of the 
original integrand, which is |.M 7 | 2 in our case, so that the phase-space integration of the 
difference can be performed numerically, even over the singular regions of the original 
integrand. In this difference, M.^ can be evaluated without regulators for soft or collinear 
singularities, i.e. we can make use of the results of the previous section. The auxiliary 
function has to be simple enough so that it can be integrated over the singular regions 
analytically, when the subtracted contribution is added again. This part contains the 
singular contributions and requires regulators, i.e. photon and electron masses have to be 
reintroduced there. Specifically, we have applied the dipole subtraction formalism, which 
is a process-independent approach that was first proposed [HU] within QCD for massless 
unpolarized partons and subsequently generalized to photon radiation of massive polarized 
fermions [l41j. We only need the limit of small fermion masses, which was worked out in 
Refs. an< ^ m wn i° n the application of the method is relatively simple. In order 

to keep the description of the method transparent, we describe only the basic structure 
of the individual terms explicitly and refer to Refs. [ El for details. 

In the dipole subtraction formalism the subtraction function is constructed from con- 
tributions that are labelled by ordered pairs ab of charged fermions, so-called "dipoles" . 
The fermions a and b are called emitter and spectator, respectively, since by construction 
only the kinematics of the emitter a leads to collinear singularities. Since we only have 
charged particles in the initial state, the subtraction function receives two contributions, 
which both have emitter and spectator in the initial state, 



\M a suh (p 1 ,P2,k i ,k)\ 2 



E \M a suh , ab (pi,p 2 ,k l ,k)\ 2 , 



a, 6=1, 2 





2 




with the dipole function 



9ab (Pa,Pb,k) 



i r 2 



(2.29) 



(p a k)x a b L 1 x a b 



1 - %ab 



where 



PaPb -Pak-pbk 



(2.30) 



•^ab 



PaPb 
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The modified momenta in Eq. ()2.28|) depend on ab and are denned as follows. While the 
spectator momentum p b is kept fixed, the emitter momentum p a , is rescaled by x ab , 



Pa ~ X abPai 



~u u 
Pb = Pb, 



and all other momenta L are transformed with a Lorentz transformation, 



kj — A^ykj , 



where 



a/j, ^ a _ (fab + PabY(Pgb + Pgb)v 2P^ b P ab ^ 



Pab Pab Pab 



p2 
r ab 



Pab = Pa+Pb- k = J2 k ji P ab = X abPa + 



(2.31) 

(2.32) 



;2.33) 



The subtracted contribution can be integrated over the (singular) photonic degrees of 
freedom up to a remaining convolution over x = x ab . In this integration the regulators 
m 7 and m e must be retained, and the soft and collinear singularities appear as logarithms 
in these mass regulators, 



* suh ( Pl ,P2,P-,P+) = ^ f'dx J2 Si snh \s 

2tt Jo £t ± 



J da {xPi, P 2,TP-,P+ 
+ / da (p 1 ,xp 2 ,P-,TP + 
+ ^ E Gi T Ub) (s) J da (s,TP-,P + ) + J da (s,P_,TP, 



(2.34) 



with the universal functions 
/l + x 2 \ 



Q ( : uh) (s,x) 



G { * nb \s) 



1 — x 

,2 



In 



>(sub) 



(s,x) = (1 - x)+, 



(2.35) 



The two cases r = +/— correspond to collinear photon emission without/with a spin-flip 
of the electron or positron. The (...)+ prescription is defined as usual, 



dx[f(x))g(x) = [ dxf(x)[g(x)-g(l)], 



(2.36) 

and P± are the degrees of polarization of the e beams. 

In summary, the phase-space integral ()2.26|) in the dipole subtraction approach reads 



a. 



^/ d$ 7 E i(l + 2P_(7)(l-2P 4 



(J 



E i^; A i 2 -i-Mr v 



sub I 



A=±l 



+ a suh . (2.37) 
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(ii) The phase-space-slicing approach 

The idea of the phase-space-slicing method is to divide the bremsstrahlung phase space 
into singular and non-singular regions, then to evaluate the singular regions analytically 
and to perform an explicit cancellation of the arising soft and collinear singularities against 
their counterparts in the virtual corrections. The finite remainder can be evaluated by 
using the usual Monte Carlo techniques. For the actual implementation of this well-known 
procedure we closely follow the approaches of Ref. [EH]- We divide the five-particle phase 
space into soft and collinear regions by introducing the cut-off parameters AE and Ad, 
respectively. We decompose the real corrections as 



dcx 7 = dcr S oft + ddcoii + da. 



7, finite- 



(2.38) 



Here d<7 so ft describes the contribution of the soft photons, i.e. of photons with energies 
ko < AE in the CM frame, and dcr co u describes real photon radiation outside the soft- 
photon region (ko > AE) but collinear to the e beams. The collinear region consists of 
the two disjoint parts < 9 7 < AO and n — AO < 1 < n, where # 7 is the polar angle of 
the emitted photon in the CM frame. The remaining part, which is free of singularities, 
is denoted by definite- 

In the soft and collinear regions, the squared matrix element |.M 7 | 2 factorizes into 
the leading-order squared matrix element |A^o| 2 and a soft or collinear factor. Also the 
four-particle phase space factorizes into a three-particle and a soft or collinear part, so 
that the integration over the singular part of the photon phase space can be performed 
analytically. 

In the soft-photon region, we apply the soft-photon approximation to lA^^,] 2 , i.e. the 
photon four-momentum k is omitted everywhere but in the IR-singular propagators. In 
this region d<7 7 can be written as [ I27| 133] 



da. 



soft 



-dcr 



a 



47T 2 



d 3 k 



A* 
Pi 



a 
P2 



k <AE 
k| 2 =fc2_ ra 



h \Pik p 2 k t 



(2.39) 



The explicit expression for the soft-photon integral can be found in Refs. [|23l2Zj- F° r 
our purpose it is sufficient to keep the electron mass only as regulator for the collinear 
singularities. In this limit we obtain 



soft 



-dcr - J 2 lis 

7T 



/2AE\ 

V m 7 / 




In the collinear region, we consider an incoming e T with momentum pi being split into 
a collinear photon and an e T with the resulting momentum xpi after photon radiation. 
In the asymptotic limit, |.M 7 | 2 factorizes into the leading-order squared matrix element 
|Alo| 2 and a collinear factor describing collinear initial-state radiation, as long as A# is 
sufficiently small. In the collinear region also the four-particle phase space factorizes into 
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a three-particle phase space and a collinear part, so that the cross section for hard photon 
radiation (k > AE) in the collinear region reads 



<rcoa(pi,P2,P-,P+) 



a 



2tt Jo 



dx J2 Q- 

T=± 



(coll). 



S,X) 



d(To(xp 1 ,P2,TP-,P+ 

+ J da (puxp 2 ,P-,TP + 



where 



,(coll), 



s,x) 



1 + x 2 
1 — X 



In 



sA6 2 
4m 2 



'(coll). 



s,x) 



X. 



(2.41) 



(2.42) 



Subtracting the soft and collinear cross sections ()2.4())1 and 1)2.41)1 from the cross section 
of the bremsstrahlung process yields the finite cross section dcT 7i fi n i te . As usual in the 
phase-space-slicing approach, this subtraction is done in practice by imposing cuts on the 
bremsstrahlung phase space, i.e. a photon-energy cut, k > AE, and a cut on the angles 
between the photon and the beams, AO < 6 1 < n — A9. 

2.3.3 Initial-state radiation in 0{a) and beyond 

In 0(a) the effect of ISR is entirely contained in the radiative corrections described 
above. However, the emission of photons collinear to the incoming electrons or positrons 
leads to corrections that are enhanced by large logarithms. In order to achieve an accuracy 
at the few 0.1% level, the corresponding higher-order contributions, i.e. contributions 
beyond 0(a), must be taken into account. This can be done in the structure-function 
method [H^JHni- According to the mass-factorization theorem, the leading-logarithmic 
(LL) initial-state QED corrections can be written as a convolution of the lowest-order 
cross section with structure functions, and the corresponding differential cross section 
reads 



J dcrisR,LL = J dxi J dx 2 T^(x 1 ,Q 2 )T^(x 2 ,Q 2 ) J da {x 1 p 1 ,x 2 p2)- 



(2.43) 



Here x\ and x 2 denote the fractions of the longitudinal momentum carried by the incoming 
electron and positron momenta just before the hard scattering process occurs. This means 
that the incoming momenta p± before emission of the collinear photon are rescaled by 
xi t 2, and the CM frame of the hard scattering process with the lowest-order cross section 
do~o(xipi, x 2 p 2 ) is boosted along the beam axis. The LL structure function including 
0(a 3 ) terms is given by 



exp (-§&7to + |&) /? 

r(i + |/C 



X) 



S 2 ( 1 4- Sr 2 

§{T^ ln(x)+4(1 + x)ln(1 



X 



[1 + x) 



5 + x 
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PI 

384 



'l + x) 



6Li 2 (x) + 121n 2 (l - x) - 3tt 5 



1 + 8a; + 3a; 2 ) ln(x) + 6(x + 5)(1 - x) ln(l - x) 



+ 12(1 + x 2 )ln(x) ln(l 

+ i(39 - 24x- 15x 2 ) 
4 



x) 



-(1 + 7a; 2 ) In 2 (a;) 



with 



ft 



7T 



and the leading logarithm 



In 



Q 2 



(2.44) 
(2.45) 
(2.46) 



Note that the scale Q 2 is not fixed within LL approximation, but has to be set to a typical 
scale of the underlying process; for the numerics we use Q 2 = s. In (|2.44|) j E is the Euler 
constant and T(y) the gamma function, which should not be confused with the structure 
functions. Note that some non-leading terms are incorporated, taking into account the 
fact that the residue of the soft-photon pole is proportional to L — 1 rather than L for 
the initial-state photon radiation. 

We add the cross section (|2.43|) to the one-loop result and subtract the lowest-order 
and one-loop contributions do"isR,LL,i already contained within this formula, 



da- 



ISR,LL,1 



dx x dx 2 [5(1 - Xl )5(l - x 2 ) + r£'>i, Q 2 )6(l - x 2 ) 
+ 5{l-x 1 )T^ 1 {x 2 ,Q 2 )} f da {x lPl ,x 2 p 2 ), 



(2.47) 

in order to avoid double counting. The one-loop contribution to the structure function 
reads 

/3 e fl + x 2 \ 



iLL,l, 



x,Q 2 ) 



X 



— hm 

4 



5(1- x) Q + 21ne) +6(1 



x 



l + x 2 



1 — X 



'2.48) 



Note that the uncertainty that is connected with the choice of Q 2 enters only in 0(a 2 ), 
if all 0(a) corrections, including constant terms, are taken into account. 

Finally, we complete the definition of the IBA by dressing dcr"B^" photomc of (|2.2(Jj) with 
the ISR structure functions, 

/ cI^ba = f Q daa £ da; 2 T^(x u Q 2 )T^(x 2 , Q 2 ) J da^ photo ™ \x lPl , x 2 p 2 ) . (2.49) 

2.3.4 Monte Carlo integration 

The phase-space integration is performed with Monte Carlo techniques in both com- 
puter codes. The first code employs a multi-channel Monte Carlo generator similar to 
the one implemented in RacoonWW [021 EZj and Lusifer [Ej, the second one uses the 
adaptive multi-dimensional integration program VEGAS [0S]. 
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3 Numerical results 

3.1 Input parameters 

For the numerical evaluation we use the following set of SM parameters [ 149] , 





= 1.16639 x lO^GeV -2 , 


or(0) 


= 1/137.03599976, 






M w 


= 80.423 GeV, 


M Z LEP 


= 91.1876 GeV, 


pLEP 
1 Z 


= 2.4952 GeV, 


m e 


= 0.510998902 MeV, 


«V 


= 105.658357 MeV, 


m T 


= 1.77699 GeV 


m u 


= 66 MeV, 


m c 


= 1.2 GeV, 


m t 


= 174.3 GeV, 


m d 


= 66 MeV, 


m s 


= 150 MeV, 


m h 


= 4.3 GeV, 



(3.1) 

which coincides with the one used in Ref. [|2T]. Since we employ the on-shell renormal- 
ization scheme, the weak mixing angle is fixed by (|2.11|) . 

As discussed in Section 12.2.41 and if not stated otherwise, we evaluate amplitudes in 
the so-called G M -scheme, i.e. we derive the electromagnetic coupling a = e 2 / (4ir) from the 
Fermi constant according to (|2.18jl . This procedure, in particular, absorbs all sizeable 
mass effects of light fermions other than electrons in the coupling a^, and the results are 
practically independent of the masses of the light quarks. The masses of the light quarks 
are adjusted to reproduce the hadronic contribution to the photonic vacuum polarization 
of Ref. [ 15 j . In the relative radiative corrections, we use a(0) as coupling parameter, 
which is the correct effective coupling for real photon emission. We do not calculate the 
W-boson mass from but use its experimental value as input. 

As explained in Section 12.2.31 we employ a fixed width in the resonant Z-boson prop- 
agator in contrast to the approach used at LEP to fit the Z resonance, where a running 
width is taken. Therefore, we have to convert the "on-shell" values of M% EP and T^ EP , re- 
sulting from LEP, to the "pole values" denoted by M z and r z in this paper. The relation 
of the two sets of values is given by [ 151] 

M z = M^ EP /-/l + (r£ EP /Af£ EP ) 2 = 91.1535 GeV, 

r z = rL EP / yfl + (rLEP/ M L EP )2 = 2.4943 GeV, (3.2) 

i.e. the difference is formally of two-loop order and numerically hardly visible in the results 
presented below. 

3.2 Results on total cross sections 

In Ref. [ E] we have already discussed numerical results for the total cross sections of 
the process e + e~ — > uuB. and the corresponding radiative corrections relative to the pure 
lowest-order prediction. Particular attention has been paid to the individual contribu- 
tions of the various sources of corrections, such as the contributions of the closed fermion 
loops, of the ISR effects in 0{a) and beyond, and of the remaining bosonic corrections. 
In this discussion the ZH-production and WW-fusion channels have been considered sep- 
arately, thereby revealing characteristic differences of the two channels. In the following 
we continue the discussion of Ref. [E]. We always sum over all three neutrino species, 
i.e. over the processes e + e~ — > i/ e i> e H, u^p^H, and v T v T H. Besides the full cross section, 
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Figure 11: Improved Born approximation (IB A) for the total cross sections and radiative 
corrections relative to the IBA for Mh = 150 GeV 



denoted "total" in the plots, in some cases we also give the cross section resulting from 
the ZH-production and WW-fusion channels separately, which are referred to as "ZH" 
and "WW" contributions, respectively. In the ZH-production channel we sum over the 
relevant contributions of all z/z/H final states, which is equivalent to multiplying the cross 
sections for e + e~ — > i/^P^H by a factor 3. 

Figure HP shows the total cross sections for the ZH and WW channels, as well as their 
incoherent ("ZH+WW") and coherent ("total") sums, in improved Born approximation 
(IBA) together with the radiative corrections relative to the IBA as functions of the centre- 
of-mass (CM) energy yfs for the fixed Higgs-boson mass Mh = 150 GeV. Analogous 
results have been given in Ref. [|2T] for the pure lowest-order cross sections cr tree and the 
corresponding relative corrections (see Figs. 1 and 2 there). While the absolute cross 
sections oiba and <7 tre e look qualitatively similar, the relative corrections normalized to 
the ctiba are systematically smaller as compared to a normalization to the pure lowest- 
order cross section. This is mainly due to the dominance of the ISR corrections which 
are properly taken into account by the IBA. In the WW channel the IBA describes 
the corrected cross section within 1% for CM energies up to ~ 500 GeV and remains 
still good within 2-3% up to ^/s < 1 TeV. At high CM energies the IBA misses non- 
universal bosonic corrections, the size of which grows with energy (compare Fig. 3 of 
Ref. [IS]). In the ZH channel the IBA deviates from the corrected cross section by <• 3% 
for y/s < 500 GeV, but even this reasonably good approximation results from accidental 
compensations between fermionic and bosonic non-ISR corrections, which are both of the 
order of 5-10% but of opposite sign (see Fig. 3 of Ref. [EI])- Above a CM energy of 
500 GeV, where ZH production is suppressed, the IBA becomes even worse, since the 
dominating Sudakov logarithms, such as aln 2 (s/M^) are missing in the IBA. 
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In Figure ^] the cross sections and their corrections are shown as functions of the 
Higgs-boson mass Mh for the typical LC energy of y/s = 500 GeV. For this CM en- 
ergy the contribution of WW fusion is much larger than the one of ZH production for 
M H < 300 GeV. For M H > y/s - M z « 409 GeV the ZH-production channel is strongly 
suppressed since the Z boson cannot become resonant there. The total cross section falls 
below 1 fb in this domain. For Mh < 350 GeV the relative corrections to the total cross 
section and to the WW-fusion channel are of the order —10% to —20% if normalized to 
Otreej but only at the level of 2% if normalized to ctiba- For the ZH-production channel 
the reduction of the relative corrections with respect to ct IB a is similar. The spikes at 
Mh = 2Mw, 2Mz, 2m t result from thresholds. These singularities are avoided if the finite 
widths of the unstable particles are taken into account appropriately (see, for instance, 
Ref. [E5). 

3.3 Results on distributions 

The distributions in the Higgs-boson energy En (defined in the CM frame) are depicted 
in Figure [T3] for y/s = 500 GeV and M H = 150 GeV. The broad continuous distribution 
for En < 250 GeV is almost entirely due to WW fusion, while the resonance peak at 
-Eh ~ 265 GeV corresponds to the Higgs-boson energy in the 2^2 kinematics of the ZH- 
production process with an on-shell Z boson. This explains the behaviour of the radiative 
corrections, which are dominated by ISR. Since the WW-fusion cross section rises with 
energy, the continuous part receives negative corrections. For the WW-fusion channel 
these are particularly large for large En near the phase space boundary. On the other 
hand, the Z peak is reduced by ISR and produces a radiative tail for En values below the 
peak, since ISR effectively reduces the scattering energy of the subsequent ZH-production 
process. This radiative tail leads to corrections of up to 220%, as can be seen in the 
inset of the lower left plot in Figure ED The relative corrections are reduced drastically if 
normalized to the IBA, which takes care of the large ISR effects. For En ^5 250 GeV the 
corrections to the total improved-Born cross section vary only weakly with En- 

Figure El illustrates the distribution in the cosine of the Higgs-boson production angle 
9n (defined in the CM frame) for y/s = 500 GeV and Mh = 150 GeV. The peak behaviour 
in the very forward and backward directions is due to the dominant WW contribution, 
while ZH production follows a shape roughly proportional to 1 — cos 2 6 i h. The relative 
corrections to the WW contribution depend only weakly on cos^h and reflect the same 
reduction after normalization to the IBA as the integrated cross sections. The relative 
corrections to the ZH contribution become large in the forward and backward directions, 
where the corresponding lowest-order cross section is small. 

Finally, we show the spectra of the photon energy E~ and the photon polar angle 
cos# 7 (both defined in the CM frame) of the radiative process e + e~ — > uuB. + 7 for 
y/s = 500 GeV and Mh = 150 GeV. In order to make the photon visible, we impose 
angular and energy cuts of 

0(7, beam) > 1° £ 7 > 0.1 GeV. (3.3) 

For small photon energies the infrared pole, which is the same for the WW and ZH 
channels, dominates the spectrum. Thus, WW fusion dominates for small E 7 , since 
the lowest-order cross section is much larger for WW fusion than for ZH production. 
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Figure 12: Cross sections and relative corrections as function of the Higgs-boson mass for 
yfi = 500 GeV 



24 



da 
dE H 

10 



fb 
GeV 



1 - 



0.1 - 



i r 

total 



full 

tree 



y/s = 500 GeV 
M H = 150 GeV 

J I I L 




da 
dE H 



fb 
GeV 



10 



|r full 

tree 



0.1 r 



0.01 r 



"i 1 1 1 1 r 

ZH production 




y/s = 500 GeV 
M H = 150 GeV 



0.001 

140 160 180 200 220 240 260 280 140 160 180 200 220 240 260 280 



En [GeV] 



En [GeV] 




-50 



yfi = 500 GeV 

M H = 150 GeV 
J I I L 



_L 



140 160 180 200 220 240 260 280 
En [GeV] 



da 

d<7lBA 

3 



2 - 
1 - 



-2 
-3 - 
-4 - 



"i r~ 

total 
ZH 
WW 




M H = 150 GeV 
y/s = 500 GeV 



_L 



_L 



_L 



140 160 180 200 220 240 260 280 
En [GeV] 



Figure 13: Distribution in the Higgs-boson energy En and corresponding relative radiative 
corrections for y/s = 500 GeV and Mn = 150 GeV 
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Figure 14: Distribution in the cosine of the Higgs-boson production angle #h and corre- 
sponding relative radiative corrections for y/s = 500 GeV and Mh = 150 GeV 
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Figure 15: Distribution in the photon energy E 1 and in the photon polar angle cos # 7 in 
the radiative process e + e~ — > uuB. + 7 for y/s = 500 GeV and Mr = 150 GeV 

For larger E^, small-angle photon emission is dominant as some kind of remnant of the 
collinear pole for forward and backward emission. This means that larger E^ reflect the 
lowest-order cross sections at smaller energies. Since the ZH-production cross section rises 
with decreasing CM energies for ^/s < 500 GeV, but the WW-fusion cross section falls off 
steeply, ZH production takes over the leading role for E y > HQ GeV. 

3.4 Comparison to related work 

Adapting the input parameters and the parametrization of the lowest-order matrix 
element to those used by Belanger et al. [US], we reproduced the numbers for the total 
cross section given in Table 2 of the first paper of Ref. [HB|- Note that we switch off the 
ISR beyond 0(a) in this comparison. In Tabled we list for each Higgs-boson mass and 
the corresponding calculated W-boson mass the results of Ref. [HH] 1 together with our 
results. The numbers in parenthesis indicate the errors in the last digits. The agreement is 
good; the relative differences are with one exception below 10~ 4 for the total lowest-order 
cross section and below 0.3% for the corrected cross section. The corrections relative to 
the lowest-order cross section agree within 0.2%. This is of the order of the statistical 
error of Ref. [I18j. which is about 0.1%. Note that Belanger et al. use a(0) to parametrize 
the lowest-order cross section. As a consequence their relative corrections are shifted by 
3Ar « +9% compared to those in the G^-scheme. 

We have also reproduced the cos^h and E^ distributions in Figures 1 and 2 of the 
first paper of Ref. [HE|- We found agreement within the accuracy of these figures. 

1 According to F. Boudjema, the numbers for the lowest-order cross section in Table 2 of Ref. [ I18j 
have integration errors of the order of 0.2%. Table ^ contains updated numbers obtained with increased 
statistics. 
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M H [GcV] 


M w [GeV] 


0"tree [fb] 


a[fb] 


c/0"tree ~ 1 


%] 


150 


80.3767 


61.074(7) 


60.99(7) 


-0.2 


Ref. [HE] 






61.076(5) 


60.80(2) 


-0.44(3) 


this work 


200 


80.3571 


37.294(4) 


37.16(4) 


-0.4 


Ref. [HE] 






37.293(3) 


37.09(2) 


-0.56(4) 


this work 


250 


80.3411 


21.135(2) 


20.63(2) 


-2.5 


Ref. [HE] 






21.134(1) 


20.60(1) 


-2.53(3) 


this work 


300 


80.3275 


10.758(1) 


10.30(1) 


-4.2 


Ref. [HE] 






10.7552(7) 


10.282(4) 


-4.40(3) 


this work 


350 


80.3158 


4.6079(5) 


4.184(4) 


-9.1 


Ref. [HE] 






4.6077(2) 


4.181(1) 


-9.27(3) 


this work 



Table 1: Total cross section in lowest order and including the full 0(a) corrections and the 
relative corrections for = 500 GeV and various Higgs masses for the input parameter 
scheme of Ref. [ HE] 1 

When considering only fermion-loop corrections, we find agreement with the calcula- 
tions of Refs. [El EH], once the appropriate renormalization and input-parameter schemes 
are adopted. For more details on this comparison we refer to Ref. [ E] • 

4 Summary 

We have presented a calculation of the complete electroweak 0(a) radiative corrections 
to the single Higgs-boson production process e + e~ — > vvH in the electroweak Standard 
Model. For e + e~ — > z/ e z/ e H, where ZH production and W-boson fusion contribute, both 
production channels are added coherently. Two methods for the treatment of the finite Z- 
boson width have been introduced. We have taken special care to treat the contributions 
from 5-point tensor integrals in a numerically stable way. The complete single-hard- 
photon matrix elements have been taken into account, where soft and collinear singulari- 
ties are treated both in the subtraction and the phase-space slicing methods. Higher-order 
ISR corrections have been included in the structure-function approach. The phase-space 
integration is performed with Monte-Carlo techniques. 

We find that the electroweak corrections are of the order of —10% and are dominated 
by ISR corrections if the lowest-order matrix element is parametrized with the Fermi 
constant G^. The non-ISR corrections are at the level of a few per cent in the G^-scheme, 
but are of the order of 10% in other schemes. At high energies, where the WW-fusion 
channel dominates, the electroweak corrections depend only weakly on the energy and the 
production angle of the Higgs-boson. 

Although the sum of ISR and non-ISR is accidentally small in the a(0)-scheme, the G^- 
scheme is nevertheless preferable for the following reasons. In contrast to the a(0) scheme, 
it does not suffer from uncertainties arising from the hadronic vacuum polarization at low 
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energies. Since the cancellation between ISR and non-ISR corrections in the a(0) scheme 
at one loop is accidental, it cannot be expected that it still holds in higher orders. One 
the other hand, the G M -scheme resums the leading universal corrections associated with 
the running of the electromagnetic coupling and the universal corrections proportional to 
the square of the top mass. 

We have shown that the corrections to the WW-fusion channel can be described by 
a simple improved Born approximation within an accuracy of typically 1% (3%) for CM 
energies below 500 GeV (ITeV). For the ZH-production channel the improved Born 
approximation, which is simply based on ISR, approximates the corrected cross section 
within 3% up to CM energies of about 500 GeV, but becomes worse at higher energies. In 
summary, the approximation can be used to include radiative corrections at the qualitative 
level, but a precision analysis will require the inclusion of the full 0(a) correction. 

Appendix 

Standard matrix elements 

The four-dimensionality of space-time implies that the SME 7Vl,f H ' CT and jCi^ w intro- 
duced in Section |2~T1 are not all independent; there are linear relations among them. 

A simple way to derive the relations with real coefficients is provided by the following 
trick. In four dimensions the metric tensor can be decomposed in terms of four indepen- 
dent orthonormal four- vectors ni, 

^ = ?^-£<nf, (Al) 
i=i 

where ■ ni — gu- Two convenient choices (j = 1, 2) for the vectors n\ are given by 

V s V s 

n% = e^n^m^s = - 2 e a ^ s p ltP p 2 , y k jtS , (A.2) 

y/Stijtzj 

with e a P~< 6 ( e oi23 _ denoting the totally antisymmetric tensor. Inserting this decom- 
position for both j = 1,2 in all contractions between different Dirac chains according 
to rj ,!r r i "'' Q = T^ ,(T g a ^T'g u , etc., and subsequently using the Dirac equation and the 
Chisholm identity 

ie^Vyg = 7 a 7 V - c, Q V + - g^^ a , (A.3) 

reduces all ZH SME Aff H,<T to four and all WW SME M™ w to two SME. 

Further relations with complex coefficients result from the direct application of the 
Chisholm identity (|A.3J) to structures like r^ pifca and subsequently using the decomposi- 
tion 

g p° = £ (Z' 1 ) 2p p iP °, Zij = 2p iPj , p 3 = h, p A = k 2 , (A.4) 
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to separate all contractions between e tensors and Dirac chains via Y u s u t al3lS = T u p v g pa 't a ^ a . 
For j^, in particular, this leads to 

T uv e k iP1 k 2 6 = T vu gP a k lPlk2 ^ = _ 2J( T T uu ( Z -l\ + T vv ( z -l\ 1 (A. 5) 



Pi V 7 12 P2 V J 22. 

with 

X = e p ^ k > = e^v^v^k^k^. (A.6) 

In the inverse matrix (Z _1 ), the determinant det(Z) occurs, which can be identified with 
det(Z) = -16X 2 . 

Altogether, the linear relations reduce the set of ZH SME to two and the set of WW 
SME to one SME. Explicitly, the relations read 

^zh,± = r zH,±^ ZH ,+ Mj m = rf m Mf w , % = 2,... 13, (A.7) 



with 



ZH,+ 

r 2 


= s, 


ZH,- 


= o, 


.WW 
'2 


— —til, 


ZH,+ 


= 0, 


ZH,- 


= S12, 


„ww 

'3 


= o, 


ZH,+ 


= o, 


ZH,- 


1 /~~f* 
— ) 


„ww 

'4 


— 2 o, 


ZH,+ 


- 1 A 

4t 22 ' 


ZH ,- 


_ til 
2 ' 


„ww 

'5 


s 

~ 2' 


ZH,+ 


_ *21 
2 ' 


ZH- 


- 1 A* 

4*12 ' 


„ww 

'6 


1 rj* 


ZH,-f 


_ *12 
2 ' 


ZH ,- 


- 

4t 2 l ' 


.WW 
'7 


i R 

4*21 ' 


ZH,+ 


_ 1 A* 


ZH ,- 


*22 


-WW 
'8 


_ «12 




2 ' 


2 ' 


ZH,+ 

r 9 


= ta^' 


ZH- 


= -4tn, 


-WW 
'9 


= 4s, 


ZH,+ 


= -2t 2 i, 


ZH- 


112 ' 


.WW 
'10 


= — ~i^~B* , 

*12 ' 


ZH,+ 

r n 


= -2ti 2 , 


ZH,- 

r n 


= -4*12, 


r WW 

' 11 


= "4*12, 


ZH,+ 
^12 


= 7^ A* 

in 


ZH ,- 
»"12 


= -4*22, 


.WW 
'12 


= 4si 2 , 


ZH,+ 
^13 


= 4, 


ZH,- 
^13 


= 16, 


.WW 
'13 


= 16 



A = SS 12 - *12*21 - *11*22 + 4iX, 
-B = SS12 + *12*21 - *11*22 + 4LY, 

C = ss 12 - *i2*2i + *n*22 + 4iX (A.9) 

Finally, there are Fierz identities relating the ZH and WW spinor chains, supplement- 
ing the above list of relations by 

Mf w = -Mi 11 '", (A.10) 

so that all SME can be expressed in terms of A4 1 ' .In terms of Weyl-van der Waerden 
spinor products, which have been defined in Section 12.3.11 fsee also Ref. [EH]), these SME 
read 

= 2(p 2 k 2 )*(piki), MT'- = 2( Pl k 2 )*(p 2 k 1 ). (All) 
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Since the lowest-order matrix elements for right- and left-handed electrons are pro- 

A ZH -4- A ZH 

portional to M 1 ' and Mi ' , respectively, the whole virtual one- loop contributions to 
the squared matrix element are of the form 

Re {Ml {Ml)*} = [fi{s, t y , Sij ) + fi(s, t fi , Sij )X] \mT'°\\ (A.12) 

where are functions of scalar products of the external momenta. Note that the f? 
and |7Vfi H,<7 | 2 are invariant under the reflection of all outgoing momenta k,, on the plane 
spanned by the beam axis and the Higgs-boson momentum k3, while X changes its sign 
under this reflection. Consequently the contribution proportional to X drops out after 
integrating over the momenta of the final-state neutrinos, which are not observable. 
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